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ABSTRACT 

In  this  paper  we  have  given  a brief  review  of  the  Mandelbrot  set,  one  of  the  best  known  icons  of  fractals  which 


arises  from  the  iteration  of  the  complex  polynomial  of  the  form  z 2 + c . We  have  discussed  about  the  role  of  critical  points 


in  such  a study  with  the  help  of  Schwarzian  derivative. 
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1.  INTRODUCTION 

The  Mandelbrot  set  has  a celebrated  place  in  fractal  geometry,  a field  first  investigated  by  the  French 
mathematicians  Gaston  Julia  and  Pierre  Fatou  as  a part  of  complex  dynamics  at  the  beginning  of  the  20th  century.  Gaston 
Julia(1893-1978)  wrote  a paper  titled  "M'emoire  sur  Viteration  des  fonctions  rationelles"  (A  Note  on  the  Iteration  of 
Rational  Functions)  [10]  where  he  first  introduced  the  modern  idea  of  a Julia  set  as  a part  of  complex  dynamics.  In  this 
paper  Julia  gave  a precise  description  of  the  set  of  those  points  of  the  complex  plane  whose  orbits  under  the  iteration  of  a 
rational  function  stayed  bounded.  In  1978  Robert  W.  Brooks  and  Peter  Matelski  investigated  some  subgroups  of  Kleinian 
groups  [17]  and  as  a part  of  this  investigation  they  first  introduced  the  concept  of  what  we  now  call  Mandelbrot  set. 

Benoit  Mandelbrot  (1924-2010)  was  a Polish-born  French  mathematician,  who  spent  most  of  his  career  at  IBM's 
Thomas  J.  Watson  Research  Center  in  Yorktown  Height,  New  York.  He  was  inspired  by  Julia's  above  mentioned  paper  on 
complex  dynamics  and  used  computers  to  explore  these  works.  In  the  year  1977,  as  a result  of  his  research,  he  discovered 
one  of  the  most  famous  fractals,  which  now  bears  his  name:  the  Mandelbrot  set.  On  1st  March  1980  Mandelbrot  first 
visualized  this  set  [18].  He  studied  the  parameter  space  of  the  complex  quadratic  polynomials  in  an  article  that  appeared  in 
the  ' Annals  of  New  York  Academy  of  science'  [12]. 

The  Mathematical  study  of  the  Mandelbrot  set  actually  began  with  the  works  of  Adrien  Douady  and  John  H. 
Hubbard  [7]  who  established  many  of  its  fundamental  properties  and  named  the  set  in  honor  of  Mandelbrot.  Interest  in  the 
subject  flourished  over  and  many  other  well  known  mathematicians  began  to  study  the  Mandelbrot  set.  Heinz-Otto  Peitgen 
and  Peter  Richter  are  the  name  of  two  such  mathematicians  who  became  well  known  for  promoting  the  Mandelbrot  set 
with  computer  oriented  graphics  and  books  [15].  A good  account  of  developing  period  of  the  theory  of  complex  dynamics 
can  be  found  in  [2], [5], [6], [22].  The  authors  are  among  the  most  active  contributors  to  this  field. 

Mandelbrot  set  may  well  be  one  of  the  most  familiar  images  produced  by  the  mathematicians  and  other  related 
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scientists  of  the  20th  century.  It  challenges  the  familiar  notion  that  the  domain  of  simplicity  and  complexity  are  entirely 
different.  Because,  the  mathematical  formula  that  is  involved  in  the  construction  of  Mandelbrot  set  consists  of  simple 
operations  like  multiplication  and  addition,  still  it  produces  a shape  of  great  organic  beauty  and  complexity  with  infinite 
subtle  variations.  The  developments  arising  from  the  Mandelbrot  set  have  been  as  diverse  as  the  alluring  shapes  it 
generates. 

The  shape  of  the  Galaxies  broke  all  Euclidean  laws  of  the  man-made  world  and  deferred  from  the  properties  of 
natural  world.  If  one  identified  an  essential  structure  like  this,  Mandelbrot  claimed,  that  the  concept  of  Mandelbrot  set,  in 
general  fractal  geometry,  could  be  applied  to  understand  its  component  parts  and  make  postulations  about  what  it  will 
become  in  future.  For  instance  interested  reader  may  see  [19]  for  study  about  distribution  of  galaxies  in  an  observed 
universe.  In  today's  world  of  wireless  communication  many  wireless  devices  use  fractal  based  compact  and  potable 
antennas  that  pick  up  the  widest  range  of  known  frequencies  [1],  [20].  Fractal  art  is  a form  of  algorithmic  art  created  by 
fractal  objects  produced  by  repeated  iterations  of  some  mathematical  rules  and  representing  the  calculated  results  as  still 
images,  animations  etc.  The  Mandelbrot  set  can  be  considered  as  a great  icon  for  fractal  art.  Graphic  design  and  image 
editing  programs  use  fractal  to  create  beautifully  complex  landscapes  and  life-like  special  effects.  Interested  readers  can  go 
through  [3], [16]  for  finding  such  applications.  Fractal  statistical  analyses  of  forest  can  measure  and  quantify  how  much 
carbon  dioxide  the  world  can  safely  process  [14].  Fractal  geometry  may  also  be  applied  to  the  various  fields  of  medicine 
such  as  cardiovascular  system,  neurobiology,  pathology  and  molecular  biology  [4],  [9]. 

The  Mandelbrot  set,  like  most  of  the  other  fractals,  arises  from  a simple  iterative  process.  The  process  involved 
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here  is  the  iteration  of  the  non-linear  relation  Zn+,  = Zn  + c on  the  points  of  the  complex  plane.  It  turns  out  that  the 

same  relation  was  already  studied  in  the  early  20th  century  by  France  mathematicians  Gaston  Julia  and  Pierre  Fatou  which 
lead  to  the  discovery  of  the  Julia  sets.  Fike  the  Mandelbrot  set,  the  Julia  set  also  have  a fractal  structure  and  are  generated 
by  using  the  same  iterative  process  employed  in  the  generation  of  the  Mandelbrot  set  but  with  slightly  different  initial 
conditions.  Interested  reader  may  go  through  [11].  There  is  only  one  Mandelbrot  set  and  infinitely  many  Julia  sets-  each 
point  on  the  complex  plane  acting  as  a parameter  to  the  Julia  set. 

The  Mandelbrot  set,  a very  beautiful  fractal  structure  enjoys  a special  status  as  a cultural  icon.  Also,  deep 
mathematics  underlies  the  Mandelbrot  set.  Despite  years  of  study  by  brilliant  mathematicians,  some  natural  and  simple-to- 
state  questions  still  remains  un-answered.  For  example,  though  Mandelbrot  set  was  known  to  the  mathematical  community 
since  1977  due  to  the  complex  form  of  shape  its  area  was  able  to  estimated  to  be  1.0565918849  ± 0.0000000028  by 

Thorsten  Forstemann  just  in  2012  [8].  Much  of  the  re-birth  of  interest  in  complex  dynamics  was  motivated  by  efforts  to 
understand  the  stunning  images  of  Mandelbrot  set,  which  is  the  prime  objective  of  this  paper. 

The  rest  of  the  paper  is  organized  as  follows:  in  section  2,  we  provide  a review  of  preliminary  concepts  and 
definitions.  Section  3 deals  with  role  of  critical  orbit  in  determining  the  dynamics  of  the  map  fc{z)=Z 2 + c . In  section 
4,  we  have  given  a brief  discussion  why  the  Mandelbrot  set  is  a fractal  of  typical  nature.  Finally,  in  section  5,  we  have 
given  a concluding  remark  of  our  study. 

2.  SOME  PRELIMINARIES 

In  order  to  carry  our  study,  we  first  need  to  provide  some  definitions  concerning  classical  deterministic  chaotic 
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dynamical  systems  which  are  discussed  in  this  section. 

Definition  2.1:  The  orbit  of  a number  z0  under  function  / :C  — > C where  C denote  the  extended  complex 
plane  i.e.,  C = Cu{°o}is  defined  as  the  sequence  of  points 

Zi  =fk)U2  =fW  ■ • =fW=fkJ,  (2>1) 

Here,  / n denote  the  nth  iterate  of  / , that  is  , / composed  with  itself  Tl  times.  The  point  z0  is  called  the  seed 
of  the  orbit. 

For  each  point  Z0  £ C , we  are  interested  in  the  behavior  of  the  sequence  given  in  (2.1)  and  in  particular,  what 
happens  as  n goes  to  infinity. 

Definition  2.2:  A point  Zq^C  is  called  periodic  point  of  / if  fn(z0  ) = Z0  for  some  integer  Tl  >1.  The 
smallest  fl  with  this  property  is  called  the  period  of  Z0 . Thus,  the  periodic  points  of  Z0  are  the  zeros  of  the  function 

F(z0,f)  = fn(z0)-z0. 

A periodic  point  with  period  one  is  termed  as  fixed  point  of  / i.e.,  Z0  is  a fixed  point  of  f if  f (z0  ) = Z0  . 

Definition  2.3:  Letv  e C . For  any  complex  valued  function  / :D— >D  where D (Z  C , the  attracting  basin  or 
basin  of  attraction  of  v under  the  function  / is  defined  as  the  set  Ay  ( v ) of  all  seed  values  whose  orbit  limits  to  the 
point  v ,i.e. 

Af  (v)  = {zeD:  fn(z)^v  } 

Note  that  the  point  v does  not  necessarily  have  to  lie  in  D . However,  if  v e D and  / is  continuous  then  for 
a,z0  e D with  fn{z0  )—>a  implies  a is  a fixed  point  of  / , as 

/( a )= /( Um  f " ( z0  ))=  lira  / ' n+1  ( z0 ) =a 

n — >°°  ft— >°o 

Fixed  points  play  a major  role  in  the  study  of  dynamical  systems.  So,  it  is  necessary  to  give  some  special  attention 
to  the  fixed  points  whenever  they  arise.  Bellow,  we  discuss  some  types  of  fixed  points  that  may  arise  in  dynamical  systems. 

Definition  2.4  (Attracting  Fixed  point):  Let  / be  a map  from  its  domain  Dc  C in  to  itself,  then 

(i)  a finite  fixed  point  agC  is  called  an  attracting  fixed  point  of  / if  there  exist  a neighborhood  U of  a such 
that  the  action  of  / moves  any  point  z in  U other  than  a closer  to  a , i.e.  | /(  z )— a \ < | z — a\ 

(ii)  Suppose,  / (oo)  = oo  then  °o  is  called  an  attracting  fixed  point  of  / if  there  exist  a neighborhood  U (zC  of 
°o  such  that  for  any  point  z in  D C\U  other  than  oo 
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| / ( Z ) | > | Z | , i.e.  the  action  of  / is  to  move  each  point  in  Dn  U — { 00  } closer  to 00  . 

Remark  2.5:  If  a is  a fixed  point  of  a continuous  map  / , then  there  necessarily  existing  some  neighborhood 
U c:  Ay  ( a)  . The  proof  of  this  does  not  require  that  / be  differentiable  at  a , but  as  we  are  only  interested  in  specific 

differentiable  functions,  here  we  give  a proof  only  for  the  case  when  | f'(a)  | < 1. 

Theorem  2.6:  Let  f:D—>D  where  D^C  be  such  that  f(a)=a,  ae  D and|  f'{a)\<\,  then  a is  an 
attracting  fixed  point  of  / . Furthermore,  there  exist  some  £ > 0 such  that  S£  (a)n  D c Ay  (a) . 

Proof:  Since  | f\a)  | <1  , one  can  select  some  8 > 0 such  that  | f\a)  | < 8 < l.By  definition, 

/'(„)=  UmZWzZk, 

Z~Cl 

Therefore  there  exist  £ > 0 such  that  for  any  z e D — {a}  and  for  which  | z — a | < £ , we  have 

< 8 =>  | f(z)  - a | < 8 | z — a | < £ 

This  shows  that  for  points  z g S£  (a) , the  function  / moves  z closer  to  a by  a factor  of  at  leasts.  Hence,  a 
is  an  attracting  fixed  point  of  / . If  we  further  iterate  the  map  / at  z > generating  the  orbit  of  z , each  application  of  / 
takes  the  corresponding  orbit  point  a step  closer  to  a . Hence,  by  using  induction  we  may  show  that 

|/w(z)  — a | < 8n  \z~  a \ < Sn  £ — > 0 Whenever,  z G D with  | z — a | < £ . 

Thus  Se  (a)  n D c Ay  (a) 

Remark  2.7:  From  above  it  is  clear  that  smaller  the  value  of  8 implies  faster  the  convergence  of  fn{z)  towards 
a . As  | < 8 , one  can  choose  the  value  of  8 according  as  the  value  of  |/  \a)\  In  particular,  if  / (a)  = a and 

f'(a)  = 0 , then  the  value  of  8 can  be  taken  to  be  extremely  small  leading  to  very  fast  convergence.  Hence,  in  such  a case 
the  fixed  point  a is  called  supper  attracting. 

Definition  2.7:  (Repelling  fixed  point)  let  / be  a map  from  its  domain  Dc  C in  to  itself,  then 

(i)  a finite  fixed  point  ae  C is  called  an  repelling  fixed  point  of  / if  there  exist  a neighborhood  U of  a such  that  the 

action  of  / moves  any  point  z in  U other  than  a further  from  a , i.e.  | /(  z )— a \ > | z — a\ 

(ii)  Suppose,  / (00)  = 00  then  00  is  called  an  repelling  fixed  point  of  / if  there  exist  a neighborhood  U cC  of  00  such 
that  for  any  point  z in  D C\U  other  than  00 


f{z)  -a 

f(z)~f{a) 

z - a 

z - a 
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/ ( Z ) < \Z  , i.e.  the  action  of  / is  to  move  each  point  in 


D n U — { °°  } further  from  oo  . 


Theorem  2.8:  Let  / : D — > D where  DcC  be  such  that  f(a)=a , as  D and  f\a)  > 1 , then  a is  a repelling 
fixed  point  of  / . Furthermore,  there  exist  some  £ > 0 such  that  for  all  z g D n S£  (a)— {a}  the  orbit  fn{z)  eventually 
leaves,  i.e.  there  exists  N such  that  fN{z)&  S£(a). 

This  theorem  can  be  proving  by  some  quick  modification  of  the  proof  of  the  Theorem  2.6. 

Definition  2.9:  The  multiplier  (or  eigenvalue , derivative)  X of  a rational  map  / iterated  n times,  at  the 
periodic  point  Z0  is  defined  as: 


X = \ 


r (z0). 

i 


/ ? 


r (z0) 


if  z0*°° 

if  z0=  00 


Where  /"  (z0)  is  the  first  derivative  of  fn  with  respect  to  z at  Z0  . 

Note  that,  the  multiplier  is  same  at  all  periodic  points  of  a given  orbit.  Therefore,  it  can  be  regarded  as  multiplier 
of  the  periodic  orbit. 

The  absolute  value  of  the  multiplier  is  called  the  stability  index  of  the  periodic  point.  It  is  used  to  check  the 
stability  of  periodic  points. 

Definition  2.10:  A periodic  point  Z0  is  called  attracting  periodic  point  if  | X | <1,  supper  attracting  if  | X | = 0 
and  is  repelling  if  X > 1 . It  is  called  indifferent  or  neutral  when  \ X = 1 


Definition  2.11:  A dynamical  system  / is  called  chaotic  if  the  following  three  conditions  are  fulfilled: 

• Periodic  points  of  / are  dense , 

• The  function  / is  transitive , and 

• / Depends  sensitively  on  initial  conditions. 


The  periodic  points  of  / are  called  dense  if  for  any  periodic  point  pl  of  / and  for  any  £ > 0 , however  small 
may  be,  the  open  sphere  S£  (p{ ) contains  another  periodic  point  p2  of  / . 

The  function  / is  called  transitive  if  for  any  pair  of  points  x and  y and  for  any  r > 0 there  is  a third  point 
ze  Sr  (x) , i.e.  the  open  sphere  center  at  x and  radius  r , whose  orbit  comes  within  the  sphere  Sr(y) . Further,  / depends 
sensitively  on  initial  conditions  if  there  is  a R > 0 such  that  for  any  S£  (x)  there  is  ye  S£  (x)  and  a positive  integer  k 

such  that  the  distance  between  / k (x)  and  / k (y)  is  at  least  R . 
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To  carry  out  our  study  for  the  rest  of  this  paper  we  consider  the  maps  of  the  form: 

fc{z)=Z2+C  (2.2) 


For  different  values  of  the  parameter  c E C Below  we  have  given  a brief  review  of  the  reason  for  choosing  such 


maps: 


Consider,  h (z)  = <2  z + ft  where,  (2^0 

Now,  hr1  {fc{h(z)))  = h~l{a2  z2  + 2a  (3  z + ft2  +c) 


(2.3) 


-az2  +2/3  z + 


j82-j3+c 


a 


By  choosing  appropriate  values  of  OC , /?  and  c we  can  make  this  expression  in  to  any  quadratic  function  / that 
we  please. 

Then 


h~\fc.h  = f=>  h~\fk.h  = fk  for  all  k 

This  means  that  the  sequence  of  iterates  {/*(z)}  of  a point  z under  / is  just  the  image  under  h~l  of  the 
sequence  of  iterates  { fk  (/z  (z))}  of  the  point  h(z)  under  fc . The  map  h transform  the  dynamical  picture  of  / to  that 

of  fc  . In  particular,  fk  (z)— » °°  if  and  only  if  fck{z)—>°°.  The  transformation  h is  called  a conjugacy  between  / and 
f c for  some  c . Any  quadratic  function  is  conjugate  to  f c for  some  c . So,  by  studying  the  dynamics  of  f c for  ce  C we 
effectively  study  the  dynamics  of  all  quadratic  polynomials.  Now,  we  define  the  basin  of  attraction  B j (°o),  filled  in  Julia 
set  K\fc),  Julia  set  J c of  the  map  fc . 

Definition  2.12  : The  set  K ( fc ) of  all  those  points  of  C which  do  not  converge  to  oo  under  iteration  of  the  map 
fc  is  called  th z filled  in  Julia  set  of  the  map  fc  i.e., 


K(fc)=  C-Bf  (oo) 

Clearly  K{fc)  is  the  set  of  all  those  points  of  C whose  orbits  are  bounded  under  iteration  of  the  map  fc  . 
Definition  2.9  : The  Julia  set  J c of  the  map  fc  is  the  boundary  dK(fc ) of  the  filled  in  Julia  set  K(fc ) . 

Note  that,  the  Julia  set  J c separates  the  two  sets  filled  in  Julia  set  K (fc)  and  basin  of  attraction  B y (oo)  of  infinity. 
Thus,  for  each  z0  e Jc,  there  is  an  open  sphere  Sr  (z0  ) with  centre  at  Z0and  radius  r > 0,  containing  a point 

U E Sr{z0)  such  that  iterates  of  U under  fc  converge  to  infinity  as  well  as  another  point  V E Sr  ( Z0  ) such  that 
iterates  of  V under  fc  do  not  converge  to  infinity. 
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Theorem  2.1  : The  filled  in  Julia  set  K (fc ) is  contained  inside  the  closed  disc  of  radious  max{|c|  , 2 }.  That  is 
K{fc ) e {z:\z  \ < max{|c|  , 2}} 

The  proof  of  this  theorem  is  immediate  consequence  of  the  following  two  lemmas. 

Lemma  2.2  : If  | C | < 2 , then  the  orbit  of  the  points  lie  outside  the  circle  of  radius  2 i.e.  the  set  of  points 
{ Z-  | Z | > 2},  escape  to  infinity. 

Note  that  if  | c | < 2 , then  K(fc  ) is  a subset  of  the  closed  disc  with  center  at  0 and  radius  2 , i.e.  if  | c |<2  then 
K{fc)z{z:\z\<2}. 

The  proof  of  this  lemma  can  be  found  in  [1 1]. 

Lemma  2.3  : If  | z |>|  c | >2  , then  | /c”  ( z ) | — >°°  as  which  means  that  z£  K(fc  ),  i.e.  ze  ( 00  ) . 

From  this  lemma  it  is  clear  that  if  | c | > 2 , then  K(fc  )cz  { z : 1 Z | < | c |} . 

The  proof  of  this  lemma  can  be  found  in  [5]. 

3.  ROLE  OF  CRITICAL  ORBIT  IN  DETERMINING  THE  DYNAMICS  OF  fc 

The  structure  of  the  Julia  set  is  strongly  influenced  by  the  behavior  of  the  critical  point  (see  definition  3.1)  of  fc. 
Clearly,  fc  has  a single  critical  point  at  z = 0 . The  orbit  of  the  critical  point  is  called  the  critical  orbit  of  fc  . The  critical 

point  of  fc  is  the  point  z for  which  the  pre-images  of  any  given  neighborhood  of  z under  / -1  are  not  all  distinct  and 
thus,  allowing  a well  defined  inverse  to  be  specified. 

Let  /(z)  be  an  analytic  map  on  C whose  power  series  expansion  at  z0  e C has  the  form 
f{z)  = f{z0)+ak{z-z0)k  +ak+l{z-z0)k+1  + . . where  ak*0 

In  this  case  we  called  z0  maps  to  /(z0)  with  degree  or  multiplicity  Vj(z0)— k . Now,  z0  is  called  a critical 
point  of  / if  v y (z0 ) > 1 . 

From  the  expansion  of  /(z) , 

0)  / u-i,  t u 

= ak{z-z0)  +ak+l[z-z0)  +.  . . 

z - z0 

Since,  k > 1 taking  limit  as  z — > z0  implies  f'(z0 ) = 0 . 

Definition  3.1:  A point  Z0  £ C is  called  critical  point  for  the  analytic  function  / if  f\z0  ) = 0 . 

3.2  The  Schwarzian  Derivative:  The  Schwarzian  derivative,  named  after  the  German  mathematician  Hermann 
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Schwarz,  is  a strange  operator  that  is  invariant  under  linear  fractional  transformation.  Rather  than  trying  to  motivate  its 
origin  further,  we  simply  define  it  and  try  to  explore  its  properties  to  find  out  the  role  of  critical  point  in  the  dynamics  of 
polynomial  maps  like  fc  . 

Definition  3.3  The  Schwarzian  derivative  of  a function  / is  defined  as 


Sf(x)  = 


D3f{x ) 3 f D2f(x) 


Df{x)  2 1 Df(x) 


Where,  Dn  f(x)  represent  the  nth  derivative  of  the  function  /(x),  n — 1, 2, 3 . 

Proposition  3.4:  If  / is  a polynomial  of  degree  at  least  two  such  that  all  roots  of  its  derivative  f'  are  real,  then 
Sf  (x)  < 0 . In  particular,  if  all  roots  of  / are  real,  then  Sf  (x)  < 0 . 

Proof:  Suppose,  ax , a2 , . . . , an  ; n >1  be  the  real  roots  of  f' , then  one  can  write 

/'(x)  = a(x-al)(x-a2) . . . (x-an). 

Taking  log  on  both  sides  we  get 

n 

log  f'(x)  = log  a log|x-a,.| 


• • f"U)  sr  1 

Differentiating,  — . = S 

/ U) 


“ x~ai 


f "(x)  f ”(-v) 

Differentiating  again,  — ,,  , — — y, — r 

f \x)  f V) 


V 1 

iix-aj- 


5/h)  = -E^J-^-irz3412  <0 

ttix-aj  2[f'(x)\ 


For  the  second  part,  suppose  that  the  distinct  roots  of  / are  < j32  < . . . < j3k  where  each  root  /?•  is  of 

k 

multiplicity  mi  ; 1 < i < k . Thus,  if  d is  the  degree  of  / then  = d . Applying  mean  value  theorem  to  / on  each 

i= 1 

of  the  interval  (/?•  , /3i+x ) we  can  find  a root  of  / 7 in  (/?•  ,j3i+l)  for  each  i = 1,  2,.  . k — 1 . Also,  f'  is  divisible  by 

k 

(x  — fit )/;/'  1 for  each  i . Therefore,  f'  has  at  least  (k  — 1)+^  (m-  —l)  = d—l  real  roots.  But  f'  is  of  degree  d—  1, 

i=\ 

so,  all  the  roots  of  / 7 must  be  real.  Hence  as  shown  above,  Sf  (x)  < 0 . 


Although  The  Schwarzian  derivative  does  not  interact  particularly  well  with  most  operations  on  functions  e.g., 


Impact  Factor  (JCC):  2.6305 


NAAS  Rating  3.19 


Mandelbrot  Set,  The  Mesmerizing  Fractal  With  Integer  Dimension 


9 


addition,  subtraction,  division  etc.,  functions  with  negative  Schawarzian  derivatives  have  very  interesting  dynamical 
properties  that  simplify  their  analysis.  The  main  reason  for  that  is  the  fact  that  this  property  (-ve  Schwarzian  derivative)  is 
preserved  by  composition  of  functions  and  consequently  by  iteration  of  function. 

Proposition  3.5:  Suppose  / and  g are  two  functions  and  h = f o g , then 

Sh{x)  = S/(g(x)).[Dg(x)]2  + Sg(x). 

In  particular,  if  Sf  < 0 and  Sg  < 0 then  Sh  < 0 
Proof:  From  the  usual  chain  rule: 

Dh(x)  = Df(g{x)).Dg(x) 

D2h(x)  = D2f(g(x)).[Dg{x)]2  + Dg(x).D2  g(x) 

D3h(x)  = £>3/(g(x)).[Dg(x)]3  + 3 D2  f{g{x)).D2  g{x).Dg{x)  + Df{g{x)).D3g{x). 

C,(  )_  D3f{g{x)).[Dg{x)]3  +3D2f(g{x)).D2g{x).Dg{x)+Df{g{x)).D3g(x) 

Df(g(x)).Dg(x) 


3 P2/(g(x)).[Pg(x)]2  + Pg(x).P2g(x) 

2 Df(g(x)).Dg(x) 


r>3/(g(x)) 

Df{g{x)) 


3fZ)2/(g(x)) 
2L  Df{g{x)) 


lr>g(x)]2  + 


Dg{x) 


l\ r>2g(x) 

2 \_  Dg{x ) 


= 5/(g(x)).[Dg(x)]2  + Sg(x) 

Now,  Sf  < 0 Sf(g(x))  < 0 also,  Sg  < 0 Sg(x)  < 0 which  implies  that  Sh(x)  < 0 . 


It  is  difficult  to  interpret  graphically  the  properties  that  a function  with  negative  Schwarzian  derivative  follow. 
However,  by  the  following  proposition  we  can  at  least  say  something  about  the  behavior  of  its  derivative: 

Proposition  3.6:  If  the  Schwarzian  derivative  of  the  function  / is  always  negative  then  its  derivative  cannot 
have  a positive  local  minimum  or  a negative  local  maximum. 

Proof:  Suppose  that  x0  is  a local  extremum  of  Df  . Which  implies  that  D2f(x0)  = Oand  hence  by  the 
definition  of  Schwarzian  derivative, 


Sf(x  0)  = 


D^fix  0) 

Df(x 0) 


By  our  assumption  / is  of  negative  Schwarzian  derivative,  therefore  we  must  have  either 


D/(x0 ) > 0 & D3  f(x0 ) < 0 or  Df{x0 ) < 0 & D3/(x0 ) > 0 . 
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Now  if  Df  has  a positive  local  minimum  at  x0 , then  by  definition  we  get  Df  (jt0 ) > 0 Therefore, 

D3f{x  0)<  0 . This  means  that  D2  f changes  sign  from  positive  to  negative  at  x0  . This  turns  out  that  x0  is  a maximum 
for  Df  rather  than  a minimum.  Hence,  Df  cannot  have  a positive  local  minimum. 

Similarly,  Df  cannot  have  a negative  local  maximum,  for  if  Df  has  a negative  local  maximum  at  x0 , then 

Df{x  0)  < 0 so  that  D3  f(x0 ) > 0 . This  means  that  D2  f changes  sign  from  negative  to  positive  at  x0  , meaning  that  x0 
is  a minimum  rather  than  a maximum. 

Theorem  3.7:  Suppose  that  Sf  is  always  negative.  If  x0  is  an  attracting  periodic  points  of  / , then  either  the 
immediate  basin  of  attraction  of  x0  extend  to  ± °o  , or  there  is  a critical  orbit  of  / whose  orbit  is  attracted  to  the  orbit  of 


Proof:  We  will  first  show  that  if  x0  is  a attracting  fixed  point  with  finite  immediate  attracting  basin,  then  the 
basin  contains  a critical  point. 

The  immediate  basin  of  attraction  of  a fixed  point  x0  must  be  an  open  interval,  for  otherwise  by  continuity,  we 
could  extend  the  basin  beyond  the  end  points.  Suppose,  (a  ,b)  be  the  immediate  basin  of  attraction  of  the  fixed  point  x0 
where  a and  b both  finite. 


Since  / :(a  ,b)—>  (a  ,b) , / must  preserve  the  end  points  of  (a  ,b)  and  therefore,  f(a)  and  f(b ) are  end  points 

of  f[(a,b)]. 

Case  1:  When /(a)  = f(b ),  i.e.  if  f(a)  = a =f{b)  or  f(a)  = b =f{b) 


Figure  3.1:  (a)  f(a)=a  = f{b) , (b)  f(a)=b  = f{b ) 

In  this  case  Rolle's  theorem  implies  that  there  is  a point  c£  (a  ,b)  such  that  Df(c)=  0 ,i.e.  (a,b)  contains  a 
critical  point  c of  f which  must  be  attracted  to  x0  . 

Case  2:  Suppose  f(a  ) - a and  f(b)  - b 
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Figure  3.2  ( f(a)=a  and  f(b)  = b ) 

Since  (a,b)  is  the  immediate  basin  of  attraction  of  x0  so  / can  have  no  other  fixed  point  in  (a,b)  other  than 
x0 . Clearly,  f{x)  > x on  (a,  x0 ) and  f(x)  < x on  (jc0  ,b),  since  otherwise  nearby  orbit  points  would  move  away  from 
x0  . Now,  the  mean  value  theorem  implies  that  there  is  a c e (a,x0)  for  which 

Dy(c)  = /(q)~  /Up)  = a- X o = 1 

a -x0  a - x0 

Note  that  c ^ x0  as  Df  (jc0  ) < 1 . 

Similarly,  there  is  a point  d e (x0  ,b)  such  thatDf(d)  = l . Therefore,  on  the  interval  (c ,d)  which  contains  x0 
in  its  interior,  we  have  Df(c ) = 1 , Df(x0 ) < 1 , and  Df(d)  = 1 . So,  Df  has  a local  minimum  somewhere  in  ( c ,d ). 
By  Proposition  3.6,  Df  cannot  have  a positive  local  minimum  in  (c , d)  and  so  it  must  attain  a negative  value  in  (c ,d). 
By  intermediate  value  theorem,  Df  must  take  the  value  0 , and  hence  / has  a critical  point  in  (a , b) . 

Case  3:  When  f(a)  = b and  f{b)  = a 


Suppose g=f2-  Clearly  ( a,b ) is  the  basin  of  attraction  of  attracting  fixed  point  x0  forg  . Now, 
g(a)  = f(f(a))  = f(b)  =a  and  g(b)  = f{f{b))=  f{a)=b.  Also  by  the  Proposition  3.5,  Sg  < 0 , therefore,  by  case 
2 above,  g has  a critical  point  c in  (a , b) . 

Now,  Dg{c)  = 0 Df{f(c)).Df(c)  =0^  Df{f{c))  = 0 or  Df{c)  = 0 
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Which  follows  that  one  of  c or  /(c)  is  a critical  point  of  / . As  both  of  them  lie  in  (a,b),  f has  a critical 
point  in  (a,b). 

Note  that  the  periodic  points  play  an  important  role  in  the  iteration  theory  and  hence  it  has  an  important  role  in 
non-  linear  dynamics  also.  As  the  critical  point  will  find  the  attracting  cycles  of  the  function  used  for  iteration,  an 
immediate  consequence  of  this  theorem  is  that  it  helps  us  to  explain  why  critical  points  are  used  to  plot  orbit  diagrams. 

An  immediate  application  of  this  theorem  to  the  map  fc  is  that  as  for  | z \ > max{2,  | c | },  the  orbit  of  z under 
fc  will  be  unbounded.  So,  there  is  no  infinite  basin  of  attraction.  Thus,  the  second  part  of  the  theorem  applies  since 
Sfc  < 0 . Also,  the  only  critical  point  of  fc  is  at  z = 0 , therefore,  if  there  is  an  attracting  periodic  points  for  fc , the  orbit 
of  z = 0 will  find  it. 


4.  THE  MANDELBROT  SET 

Generally,  in  the  study  of  polynomial  function  / as  a dynamical  system,  we  first  choose  a seed  value  z0  and 
then  try  to  understand  the  long  term  behavior  of  the  sequence:  z0 , z1  = f{z0 ),  z2  = f(zi  )>  • • • More  particularly,  we 
will  try  to  answer  such  question  as: 

• For  a fixed  polynomial  / what  seed  value  leads  to  bounded  sequence? 

• For  a parameterized  family  of  polynomials,  how  do  the  set  of  such  seeds  depend  on  the  parameter  ? 

Mainly  for  most  of  the  families  of  polynomial  systems,  there  is  no  obvious  picture  to  make  in  the  parameter  space, 
since  there  is  no  obvious  question  to  ask.  The  study  of  the  dynamics  of  complex  polynomials  and  rational  functions  is  a 
success  story  mainly  because  of  the  role  played  by  the  critical  points,  and  therefore,  in  this  cases  we  may  study  about: 

• What  happens  to  the  critical  points  under  iterations? 

The  Mandelbrot  set  is  an  answer  to  this  question  for  the  complex  one  parameter  family  of  quadratic  polynomials 
{fc-ce  C}. 

The  subset  of  the  parameter  plane  ( or  c-plane)  consisting  of  all  parameter  value  c for  which  the  orbit  of  the 
critical  point  z—  0 under  the  map  fc  , i.e. 


Is  bounded  is  termed  as  the  Mandelbrot  set. 

Definition  4.1:  The  Mandelbrot  set  M is  defined  as: 

M = { c e C : The  Orbit  of  0 is  bounded  under  iteration  by  fc } 

= {ce  C:  3r  > 0,  | fcn  (o)|  < r,  Vne  n} 

One  of  the  particular  interests  is  to  represent  Mandelbrot  set  graphically.  The  simplest  algorithm  for  generating  a 
representation  of  the  Mandelbrot  set  is  known  as  'Escape  Time'  algorithm  where  we  color  each  point  on  the  parameter 
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space  depending  on  where  its  attractor  lies  i.e.  whether  it  is  attracted  to  infinity  or  bounded  within  the  set. 

In  view  of  Lemma  2.3  it  is  clear  that,  if  | c | > 2 , then  | fc  (o)  | = | c | > 2 and  therefore  the  orbit  of  the  critical 
point  z = 0 necessarily  escape  to  infinity  in  this  case.  To  distinguish  these  points  we  assign  a particular  color,  say  color- 1, 
to  these  points,  i.e.  the  points  for  which  | c | > 2 . Now,  it  is  clear  that  Mandelbrot  set  constituted  by  some  points  in  the 
parameter  plane  within  the  [-  2,  2]  X [-  2, 2]  grid.  To  find  out  these  points,  we  first  choose  a maximum  number  of 
iteration,  say  N , as  the  "bailout"  limit.  Then,  for  each  c value  in  this  grid,  we  compute  the  set  of  points 
{ /cW(o)  : = — , 2V  }.  If  for  some/:  (0  < k < N),  //( 0)  > 2,  then  clearly  | /cn(z)|  — > oo  as  n —>  oo , and 

therefore,  we  stop  further  iteration  and  assign  color- 1 to  these  points.  Again  if,  | (o)  | < 2 V k = 1, 2,  • • • , N then  we 

assign  another  color  (contrasting  to  color- 1),  say  color-2  to  this  point  in  the  parameter  plane.  After  completing  the  iteration 
process  for  all  the  c - values  the  region  covered  by  color-2  is  an  approximation  of  the  Mandelbrot  set  M . 

It  is  to  be  noted  that  though  this  algorithm  always  correctly  identifies  the  points  outside  the  Mandelbrot  set,  the 
imposed  maximum  number  of  iteration  (i.e.  N ) causes  the  algorithm  to  mis-classify  points  as  being  inside  the  set  since  for 
starting  values  very  close  to  but  not  in  the  Mandelbrot  set  may  take  hundreds  or  thousands  of  iterations  to  escape.  This  mis- 
classification  of  points  near  the  boundary  of  M set  can  be  notice  by  co  mparing  two  figures  given  bellow  (Figure  4.1) 
where  Figure  (a)  and  (b)  are  generated  by  taking  N = 10  and  N = 150  respectively  and  considering  'green'  as  'color- 1' 
and  'black'  as  'color-2'. 


(a)  (b) 

Figure  4.1:  Graphical  Representation  of  Mandelbrot  Set  [(A)  N=10(B)  N=150 

As  the  Figure  4.1  (a)  is  generated  by  less  number  of  iterations,  this  figure  mis-classified  more  points  as  bounded. 

Thus,  higher  the  value  of  N , i.e.  the  number  of  iterations,  will  give  more  detail  image  of  the  Mandelbrot  set,  but  in  this 
case  the  computer  will  take  more  time  for  generating  the  image. 

The  Mandelbrot  set's  true  visual  beauty  rely  on  the  coloring  near  its  boundaries.  Developing  a strong  coloring 
algorithm  helps  display  the  beauty  of  the  set  by  providing  the  stunning  visual  aspect  of  the  set,  which  also  gives  the 
excitement  of  studying  the  set?  One  of  the  most  popular  way  of  doing  this  is  by  assigning  different  colors  to  the  points  in 
the  various  regions  such  as  inside  the  set,  boundary  of  the  set.  Also,  for  the  points  just  outside  the  boundary,  colors  are 
determined  by  the  number  of  iterations  needed  by  the  point  to  exceed  a certain  test  value  (usually  2 ).  In  the  Figure  4.2  , we 
use  blue  for  the  points  inside  the  set,  green  for  that  in  the  boundary  of  the  set  and  orange  color  for  the  points  just  outside 
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the  set.  Gradually  deeper  color  in  orange  indicates  less  number  of  iterations  needed  to  exceed  test  value  2 in  magnitude. 


Figure  4.2:  (The  Mandelbrot  Set) 

4.2  Some  Observations  on  Mandelbrot  Set 

Computer  images  of  the  Mandelbrot  set,  as  in  the  Figure  4.2  , shows  the  estimated  geometry  of  the  Mandelbrot 
set.  It  contains  a big  cardioid  shaped  region,  called  the  body  of  the  Mandelbrot  set.  This  region  is  indicated  by  B in  Figure 

1 3 

4.3  and  it  intersects  the  real  axis  at  c — — & c — . Towards  its  left,  a circular  area  H with  center  at  c — — 1 and 

4 4 

radius  — is  attached,  called  the  head  of  the  set. 

4 


Figure  4 3:  (The  Body  B and  Head  H of  the  Mandelbrot  Set) 

The  surface  of  these  two  parts  are  covered  by  some  richly  detailed  structure  of  decoration  which  makes  the  set  a 
fractal  one.  Closer  inspection  of  these  decorations  shows  that  all  of  them  are  different  in  shape.  Any  such  decoration 
directly  attached  to  the  body  is  called  a primary  bulb  or  decoration.  In  turn,  there  are  many  smaller  decorations  attached  to 
the  boundary  of  each  of  these  decoration  as  antennas.  Again,  antennas  attached  to  each  decoration  seems  to  consist  of 
several  spokes.  The  number  of  such  spokes  varies  from  decoration  to  decoration  as  clearly  visible  in  the  Figure  4.4. 
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(e)  (f) 

Figure  4.4:  (Some  Decorations  on  the  Primary  Bulb) 


The  Mandelbrot  set  is  indeed  a fractal  object,  in  the  sense  that  it  has  repeating  patterns  at  different  scales.  In  fact, 
its  boundary  is  filled  with  a halo  of  tiny  copies  of  the  entire  set,  usually  referred  to  as  a 'baby  Mandelbrot  set'.  Each  of  these 
baby  Mandelbrot  set  is  again  surrounded  by  its  own  halo  of  still  tinier  copies,  and  so  on,  smaller  and  smaller  scales  without 
end.  For  example,  Figure  4.4  shows  successive  magnifications  of  a portion  of  the  Mandelbrot  set.  Where  the  region 
indicated  by  a rectangle  in  the  Figure  4.4(a)  is  magnified  and  plotted  as  Figure  4.4(b).  Again  the  portion  of  the  Figure 
4.4(b)  covered  by  the  rectangle  is  magnified  and  plotted  as  Figure  4.4(c),  and  so  on.  Note  that,  though  these  Figures 
certainly  suggest  that  the  baby  Mandelbrot  sets  are  like  island,  well  separated  from  one  another  and  from  main  body  of  the 
set,  this  fact  is  not  true  in  reality.  Douady  and  Hubbard  have  shown  that  the  Mandelbrot  set  is  connected  [7].  That  is,  the 
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intricately  branched  decorations  are  filled  with  filaments  of  baby  Mandelbrot  set.  Though  they  are  not  visible  at  certain 
level  of  magnification,  they  all  link  back  to  the  main  body  of  the  set. 


E 

(a)  N=  42  (b)  N=  150 

(c)  N=210 

a 

4 

f 

4 |n 

(d)  N=425  (e)  N=  900  (f)  N=  9500 

Figure  4.4:  (Baby  Mandelbrot  Set) 


Though  the  Mandelbrot  set  is  a fractal  object,  its  boundary  is  so  complex  and  intricate  that  it  has  a integer 
dimension  [21].  The  observation  that  each  filament  in  the  decoration  of  the  Mandelbrot  set  is  filled  with  baby  Mandelbrot 
set  might  lead  to  the  wrong  conclusion  that  the  Mandelbrot  set  is  self  similar.  Actually,  as  Figure  4.5  suggests,  every  baby 
Mandelbrot  set  has  its  very  own  pattern  of  external  decorations,  everyone  different  from  every  other,  i.e.,  the  baby 
Mandelbrot  sets  are  not  exact  replicas  of  the  full  Mandelbrot  set.  Which  lead  us  to  conclude  that  the  Mandelbrot  set  is  not 
exact  self  similar  although  it  appears  similar  to  our  eyes.  Mandelbrot  introduced  the  term  'statistical  self  similarity'  to 
represent  such  type  of  quasi  similar  objects  [13].  Note  that,  the  points  constituting  these  (statistical  self-similar)  objects 
belongs  to  the  same  statistical  distribution. 


(a)  N = 100  (b)  N = 360 
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(c)  N=  750  (d)  N=  3400 

Figure  4.5:  (Baby  Mandelbrot  Sets  are  Not  Exact  Self  Similar) 


5.  CONCLUSIONS 

The  Mandelbrot  set  is  one  of  the  most  beautiful  examples  of  the  fascinating  world  of  fractals.  Every  little  piece  of 
it  is  loaded  with  some  beautiful  almost  self  similar  structures.  We  tried  to  explore  those  beautiful  structures.  The  study  can 
be  useful  in  teaching  and  learning  about  the  world  of  fractals. 
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